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ABSTRACT 

We study the evolution of supernova remnants in a low-metallicity medium Z/Z Q = 10~ 4 
- 10~ 2 in the early universe, using one-dimensional hydrodynamics with non-equilibrium 
chemistry. Once a post-shock layer is able to cool radiatively, a dense shell forms behind 
the shock. If this shell becomes gravitationally unstable and fragments into pieces, next- 
generation stars are expected to form from these fragments. To explore the possibility of this 
triggered star formation, we apply a linear perturbation analysis of an expanding shell to our 
results and constrain the parameter range of ambient density, explosion energy, and metallic - 
ity where fragmentation of the shell occurs. For the explosion energy of 10 51 ergs(10 52 ergs), 
the shell fragmentation occurs for ambient densities higher than > 10 2 cm -3 (10 cm -3 , re- 
spectively). This condition depends little on the metallicity in the ranges we examined. We 
find that the mode of star formation triggered occurs only in massive (> 10 8 Af Q ) haloes. 

Key words: cosmology:theory - galaxies:formation - high-redshift - ISM:supernova rem- 
nants - shock waves - stars: formation 



1 INTRODUCTION 

Understanding of star formation and its feedback effects in the 
high-redshift universe is a key to unravelling how the primeval uni- 
verse has evolved into the variety of luminous objects we observe 
today. Recent numerical studies have indicated that the first stars, 
or Pop III. 1 stars (O'Shea et al. 2008), presumably have typical 
mass of > 100 M Q (Abel, Bryan & Norman 2002; Bromm, Coppi 
& Larson 2002; Omukai & Palla 2003; Yoshida et al. 2006). Copi- 
ous amounts of radiation from such very massive stars photo-ionize 
hydrogen atoms and photo-dissociate molecular hydrogen in their 
surroundings (Kitayama et al. 2004; Whalen et al. 2004). As well 
as the radiative feedback, the first star deaths have an impact on 
subsequent star formation around them. The ultimate fate of metal- 
free stars is determined by their mass. Stellar evolution models pre- 
dict that stars with 40 - 140 Mq finally collapse into black holes 
without supernovae, while those with 10 - 40 Mq or 140 - 260 
M0 explode as Type II supernovae or pair-instability supernovae 
(PISNe), respectively (Heger & Woosley 2002; Umeda & Nomoto 
2002). 

If the first star dies without supernova, the Hll region is left 
over around a remnant black hole. At the centre of the fossil Hll 
region, H2 is quickly replenished owing to a large amount of elec- 
trons, which are the catalysts for the H2 formation reaction. The re- 
sultant fast radiative cooling allows a small portion of the gas to re- 
collapse and turn into the next-generation stars, provided that no ex- 
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ternal ultraviolet (UV) radiation effectively dissociates H2 and the 
dark halo mass is higher than 10 s Mq (Nagakura & Omukai 2005). 
In addition, deuterated hydrogen (HD) molecules abundantly form 
in such regions and further reduce the gas temperature to that of the 
CMB. The re-collapsing gas fragments with typical mass of ~ 40 
Mq (Yoshida et al. 2007), which is an order of magnitude lower 
than that of the first stars (see also, Uehara & Inutsuka 2000; Naka- 
mura & Umemura 2002). 

If the first star explodes as a supernova, the surrounding 
medium is severely affected by the blast wave. The blast wave re- 
moves most of the gas from a low-mass halo (< 10 6 Mq) even 
for the explosion energy as low as 10 50 erg (Kitayama & Yoshida 
2005; Whalen et al. 2008), because the gas density has been de- 
creased to 0.1 - 1 cm~ 3 inside the Hll region. Evacuation of the gas 
from the halo quenches subsequent star formation in the same halo 
without a substantial gas supply through merging with other haloes 
(Greif et al. 2007). For a supernova in higher-mass haloes (> 10 7 
Mq), the situations are quite different. The halo gas remains in the 
Hll region owing to a deeper dark matter potential well (Kitayama 
et al. 2004). The blast wave stagnates inside the halo and sweeps a 
considerable amount of gas into a dense shell (Whalen et al. 2008). 
Such a shell is a potential site for next-generation star formation. 
The stability analyses of linear perturbations on an expanding and 
decelerating thin shell predict that small density perturbations can 
grow to trigger gravitational instabilities, provided that the temper- 
ature in the shell becomes low enough and sufficient gas is taken 
into the shell (Elmegreen 1994; Whitworth et al. 1994). Some au- 
thors have investigated the possibility of star formation triggered 
by the first supernova explosions (Ferrara 1998; Mackey, Bromm 
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& Hemquist 2003; Bromm, Yoshida & Hernquist 2003; Salvaterra, 
Ferrara & Schneider 2004; Machida et al. 2005; Vasiliev, Vorobyov 
& Shchekinov 2008; Whalen et al. 2008). Applying a classical su- 
pernova remnant (SNR) model to a primordial gas, Salvaterra et 
al. (2004) and Machida et al. (2005) found the combination of the 
explosion energy and the ambient density where the shell fragmen- 
tation occurs, and discussed the typical mass of the stars forming in 
the shell. However, their results are somewhat different because of 
different one-zone modelling. To solve this uncertainty, more de- 
tailed calculation is needed. 



Supernovae also disperse heavy elements into the intergalac- 
tic medium (IGM) (e.g., Kitayama & Yoshida 2005). Greif et al. 
(2007) showed that the metals are preferentially transfered into 
voids, while some of them collide with neighbouring haloes and 
mix with the gas in them through hydrodynamical instabilities such 
as Kelvin-Helmholtz instability. Cen & Riquelme (2008) showed 
that the mixing occurs only at the edge of haloes of 10 6 " 7 Mq at 
z — 10. However, irradiation of ionizing photons from neighbour- 
ing haloes, which was omitted from their analysis, may enhance the 
metal enrichment of lower-mass haloes at higher redshifts. Turbu- 
lent motions excited by hierarchical mergers also causes the effi- 
cient metal mixing (Wise & Abel 2007; Grief et al. 2008). Metals 
contribute to the radiative cooling of gas, having an effect on dy- 
namical evolution of SNRs. Thornton et al. (1998) investigated this 
issue for a wide metallicity range with the intention of including the 
SN feedback in galaxy-formation simulations. They showed that 
the metal cooling modifies the redistribution of supernova energy 
to the surroundings for metallicity > O.OIZq. However, they paid 
no attention to the structure and evolution of the supernova-driven 
shell nor to the possibility of subsequent fragmentation. 



The aim of this paper is to find the condition where the 
supernova-driven shell fragments and triggers subsequent star for- 
mation in young galaxies. Towards this end, we study the evolution 
of supernova remnants in the low-metallicity interstellar medium 
by way of hydrodynamics with spherical symmetry. This paper 
is organized as follows: In Section [2] we describe the numerical 
method and input physics. The results are presented in Section [3] 
where we first describe dynamical evolution of an SNR and the 
structure of the shell, focusing on the roles of metals, and then 
discuss the possibility of the shell fragmentation using the lin- 
ear stability analysis by Elmegreen (1994). In Section [4] by using 
the above result, we constrain the mass range of haloes where the 
triggered star formation occurs. In addition, we discuss effects of 
processes not included in our analysis. Finally, we summarize our 
study in Section [5] Throughout this paper, we adopt a flat ACDM 
model with parameters S1a = 0.73, Q m = 0.27, ttb = 0.044, and 
Hubble parameter h — 0.71 (Spergel et al. 2003). 



2 THE MODEL 



2.1 Hydrodynamics and Chemistry 

We solve the following basis equations for one-dimensional hydro- 
dynamics with spherical symmetry: 
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where p is the mass density of the gas, v the velocity of the fluid ele- 
ments, p the pressure, e the internal energy per unit mass, m(r) the 
enclosed mass within the radius r, 7 the adiabatic exponent, p, the 
mean molecular weight, mn the mass of a hydrogen atom, A the 
radiative cooling rate per unit volume and time, &b the Boltzmann's 
constant, and n(T) the thermal conductivity. The thermal conduc- 
tivity is the sum of contributions from atomic diffusion at low tem- 
peratures (Parker 1953) and electronic diffusion at high tempera- 
tures (Spitzer 1962; Cowie & McKee 1977). In our calculations, a 
cooling layer behind a shock driven by a supernova becomes ther- 
mally unstable for wavelengths shorter than a so-called Field length 
(Field 1965) : X F = (kT/A) 1/2 . Thermal conduction tends to cur- 
tail the thermal instability, but the Field length is many orders of 
magnitude shorter than a characteristic length (e.g., a radius of an 
SNR). The limitation of our simulation resources makes it diffi- 
cult to resolve the Field length. To suppress unphysical oscillations 
due to the thermal instability occurring below the spatial resolu- 
tion limit, we temporarily enhance the thermal conductivity only 
for the thermally unstable cooling layer at the incipient stage of the 
shell formation (see Hosokawa & Inutsuka 2006). This prescription 
smoothes the temperature structure out on the newly defined Field 
length, but does not significantly alter the thermal evolution of the 
shell because the Field length is still shorter than the scale-length 
of temperature variation. 

The hydrodynamics are solved with a second-order Godunov 
method in the Lagrangian coordinate (van Leer 1979). We split the 
gas into fluid elements which are unequal in mass such that more 
fluid elements exist in the post-shock regions. We set the number 
of the total fluid elements to 1000 and choose a simulation box of 
length 100 - 1000 pc, depending on the explosion energy and the 
ambient density. We have confirmed that the results do not change 
when the total number of mesh points is doubled. To check the 
accuracy of our code, we have carried out a test simulation on a 
dynamical expansion of a point-source explosion without any dis- 
sipative processes and ascertained that the result is excellently con- 
sistent with the Sedov-Taylor solution. 

Along with the hydrodynamics, we solve a set of kinetic equa- 
tions for chemical species. We include 15 primordial components 
and 4 heavy elements as follows : H, H + , H~, H2, Hj, HeH + , 
He, He+, He ++ , e~ , D, D+, D", HD, HD+, C, C+, O, and 0+. 
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The kinetic equation for the i-th species is 

= -nny(i) ^ k ljy (j) + n H ^ h m y{l)y{m) 

j lm 

+ n1 I ^2k i j k y(i)y(j)y(k), (5) 

ijk 

where nu is the number density of hydrogen nuclei, y(i) the frac- 
tion of the number density for the i-th chemical species to riH, fcij 
and kim the destruction and the formation rate of y(i), respectively, 
and kijk the reaction rates for three-body reactions. We include 53 
reactions among these species on the basis of the minimum model 
suggested by Omukai et al. (2005) and listed in Table lAll The pri- 
mordial chemical networks are based on Abel et al. (1997) and 
Galli & Palla (1998). The reactions involved with carbon and oxy- 
gen are the recombination of C + (reaction [R42]) and + (reaction 
[R43]). In addition, we incorporate the charge exchange reactions 
between H and + (reaction [R44]), and between H + and O (re- 
action [R45]) because these reactions occur quickly due to a slight 
difference among their ionization potentials. We do not take into ac- 
count the charge exchange reactions between C and H + or between 
H and C + because these reactions are much less important than re- 
actions [R44] and [R45] (Glover & Jappsen 2007). Although the 
emissions from the compounds of carbon or oxygen such as OH, 
and H2O slightly contribute to the total radiative cooling at den- 
sities ranging from 10 s cm -3 to 10 10 cm -3 (see, Omukai et al. 
2005), we omit the reactions among them because they have only 
minor effects in the density range we explored. We apply the set 
of equations to each fluid element and solve them with an implicit 
scheme because of the stiffness of these equations. The adopted 
time steps are actually shorter than the chemical time. We implicity 
solve the chemical reactions to make sure the numerical stability. 
We have confirmed that our results do not change even with the 
shorter time steps. 

2.2 Radiative Processes 

The net cooling rate A in equation ((3) is the sum of the con- 
tributions from radiative cooling processes. In simulating super- 
nova explosions, we need to treat both the hot interior bubble and 
dense shell driven by the blast wave. The temperature range is 
wide, spanning a wide range from 10 K to 10 9 K. In the high- 
temperature regime of > 10 4 K, the important processes include 
thermal bremsstrahlung emission, inverse Compton scattering of 
the CMB, collisional ionization and excitation of hydrogen and he- 
lium atoms, recombination of hydrogen and helium. Their cooling 
rates are taken from Cen (1992). On the other hand, at lower tem- 
peratures, the radiative processes by molecules and metals become 
important. We take account of H2 and HD ro-vibrational cool- 
ing (Hollenbach & McKee 1979; Galli & Palla 1998). Recently, 
Glover & Abel (2008) have computed new H2 cooling rates for 
e~ - H2 and H + - H2 collision channels in addition to the rate 
we included. We do not include these updated rates for simplic- 
ity. This effect might slightly reduce temperature of the post-shock 
gas, but not affect our conclusions. We also include the meta-stable 
and fine-structure line cooling by carbon and oxygen (Hollenbach 
& McKee 1989): [Cl]A9823, A4622, and A8727 ; [Cll]A2326 ; 
[Oi]A6300, A2972, and A5577 ; [On]AA3729, 3726, and 508 /xm 
; [Cl]609.2Aim, 229.9um, and 369.0Atm ; [Ol]63.1um, 44.2/xm, 
and 145.6^im. The cooling rates are derived by solving the statisti- 
cal equilibrium between the energy levels. We omit the metal line 
cooling at temperatures higher than 10 4 K because its contribution 



is less important in the metallicity range we consider (see, Suther- 
land & Dopita 1993). 

Some metals are expected to have condensed into dust grains 
even in the early universe (Nozawa et al. 2003; Bianchi & Schnei- 
der 2007). The dust can be an efficient radiator, but it is important 
only at densities far higher than we consider here. Also, it plays 
a role in a heating source through photo-electric heating under a 
strong radiation field. This rate depends on the properties of the 
radiation field as well as dust grains, both of which are quite uncer- 
tain in the high-redshift universe. Here, we just omit the dust effect 
for simplicity. 

The CMB temperature is quite a high value of ~ 44 K at 
z ~ 15. If a gas temperature decreases below that of the CMB, the 
CMB photons heat the gas through radiative pumping. This effect 
is taken into consideration in our calculation. We do not include 
other heating such as cosmic rays or external radiation from other 
sources. 



2.3 Initial Condition 

Before the supernova explosion, the progenitor star ionizes the sur- 
rounding medium. Dynamical expansion of the Hll region has been 
investigated by several authors (Whalen et al. 2004; Kitayama et 
al. 2004; Hosokawa & Inutsuka 2006). A shock wave associated 
with the ionization front pushes a large amount of the circumstel- 
lar gas away, leveling the density off. Thus, we simply assume that 
the surrounding gas is homogeneous and fully ionized at tempera- 
ture 10 4 K. In reality, the typical temperature in Hll regions around 
Population III stars is somewhat higher than 10 4 K (Whalen et al. 
2004; Kitayama et al. 2004). However, this difference will not be so 
significant in our calculations. In an early phase of the expansion, 
ambient thermal pressure is negligible compared with ram pressure 
exerted on the shell. Temperature in the fossil Hll regions quickly 
falls below 10 4 K by the Lyman-a cooling. The initial temperature 
difference will be wiped out before the ambient thermal pressure 
acts to decelerate the shell. The ambient density depends on the 
conditions of formation sites of the progenitors. Here, we treat it as 
a free parameter and study the cases of no = 0.1, 1, 10, and 10 3 
cm -3 . We start the simulations at redshift z = 15. 

We inject the supernova energy Esn in the form of thermal 
energy and load the ejecta mass into several central meshes. Since 
the explosion energy can differ from event to event, we studied 
the three cases of £sn = 10 51 , 10 52 and 10 53 ergs, which corre- 
spond to normal core-collapse supernovae, hypernovae and PISNe, 
respectively. Note that we performed the PISN runs only for low- 
metallicity (Z = 10~ 4 Zq) because metal enrichment reduces a 
typical fragment scale of a collapsing gas, and hence stellar mass 
(e.g., Omukai et al. 2005). The result depends little on the ejecta 
mass because it is much smaller than the mass swept up later into 
the shell. Hence, we simply fix it to 10 Mq. 

For the elemental abundances in ISM, the number fractions of 
He and D nuclei to the hydrogen ones are 7.9 x 10 -2 and 2.5 x 
10~ 5 respectively. We study the cases with metallicities Z/Zq — 
10~ 4 , 10~ 3 , 10" 2 . For simplicity, we assume the metal abundance 
ratios are in proportion to the solar values. The number fractions of 
carbon and oxygen are 3.6 x 10~ 4 and 8.5 x 10~ 4 , respectively, 
for the solar metallicity (Anders & Grevesse 1989). 
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3 RESULTS 



3.1 Dynamical Expansion of a Supernova Remnant 

We first present the result for the parameter set of _Esn = 10 52 
erg, no = 1 cm -3 , and Z — 10~ 4 Zq, which we hereafter call the 
fiducial run. Fig. [T] shows the snapshots of (a) the density, (b) tem- 
perature, (c) pressure, and (d) velocity distribution as a function of 
the radius. The numbers in Fig.[T]represent the structures at time (1) 
9.0 x 10 3 , (2) 4.8 x 10 4 , (3) 1.7 x 10 5 , (4) 4.9 x 10 s , (5) 1.1 x 10 6 , 
(6) 5.0 x 10 6 , and (7) 1.0 x 10 7 yr. A contact discontinuity forms 
quite early at the boundary between the stellar remnants and the 
surrounding medium, which is consistent with the previous studies 
(Chevalier 1974; Cioffi, McKee & Bertschinger 1988; Thornton et 
al. 1998). However, this discontinuity disappears before the state 
(1) because the heat transfer causes the gas slightly outside it to 
flow into a hot interior bubble while smoothing out the boundary. 
The strong shock evolves adiabatically until t — 8.5 x 10 4 yr. The 
post-shock density is (7 + l)/(7 — 1) =4 times as high as that 
of the ambient medium. The shock front moves outwards in con- 
formity with the canonical power-law form ; oc t v , where the 
exponent r\ — 2/5 during the Sedov-Taylor (ST) phase. 

The expansion time of the shock front, i? s h/7? s h, increases 
linearly with the time elapsed during the ST phase, while the ra- 
diative cooling time in the post-shock layer, pe/A, continues to 
decrease. After the cooling time becomes shorter than the expan- 
sion time at t = 8.5 x 10 4 yr, the post-shock temperature quickly 
plunges to a low value and the gas in the post-shock layer is com- 
pressed to a thin shell. The mean pressure in the hot bubble is sev- 
eral orders of magnitude higher than that of the ambient medium 
at this moment (Fig. Q] c), and drives the outward motion of the 
shell. This is called the pressure-driven snowplough (PDS) phase. 
An analytic model of a point-source explosion exhibits r\ = 2/7 in 
the absence of radiative cooling in the hot bubble (McKee & Os- 
triker 1977). The power-law index of our results, however, deviates 
from that of the analytic model 2/7 because the interior pressure 
at the ST phase slightly increases the effective value of r\ (Cioffi 
et al. 1988; see Whalen et al. 2008, in detail). Note that the ex- 
pansion law of the bubble will change with more realistic density 
distribution in the relic HII region, e.g., dumpiness, and radial den- 
sity gradient. These effects should be separately examined in future 
studies. This phase continues as long as the mean interior pressure 
is much higher than the ambient pressure. 

Radiative cooling is inefficient inside the bubble because of 
the extremely low density, and the gas loses thermal energy through 
adiabatic expansion instead. The mean interior pressure continues 
to decrease and becomes comparable to the ambient pressure at 
t = 4.6 Myr. The motion of the shell is then driven by its momen- 
tum. This is a so-called momentum-conserving snowplough (MCS) 
phase. The thin shell model predicts that r\ — 1 /4, but it is slightly 
higher in our results because of the contribution of finite pressure in 
the bubble (Cioffi et al. 1988). The shell's expansion velocity even- 
tually falls to the sound speed of the ambient medium. In reality, 
the shell loses its identity thereafter and mixes with the ISM. This 
is the end of the lifetime of the shell. 

The results of the other runs are similar to the fiducial one. 
In all cases, the SNRs enter the MCS phase before the end of the 
shell lifetime. This is because the ambient temperature monotoni- 
cally decreases to a low value by radiative cooling. If the ambient 
matter is heated, for example, by external radiation, the MCS phase 
is unlikely to be reached. 



3.2 Thermal Evolution and Internal Structure of the Shell 

Fig.|2]shows time evolution of the internal structure of the shell. Ra- 
dial density and temperature profiles are presented against column 
density of hydrogen nuclei measured from the center. The column 
density of the shell increases as the shock front sweeps up the ambi- 
ent materials, whereas the column density of the bubble decreases 
with time, because the mean density in the bubble quickly falls by 
the expansion. The mean density in the shell gradually decreases 
over the snowplough phase. This is because the shock strength be- 
comes weaker and weaker as the bubble expands. The minimum 
temperature in the shell decreases with time. Note that it eventually 
falls below the CMB value. To see the physical processes, we plot 
in Fig.[3]the evolutionary trajectories of some typical fluid elements 
on the temperature-density diagram. These elements are taken into 
the shell in the snowplough phase. After they slightly drop from 
the initial location at (no, -Esn) = (1 cm -3 , 10 4 K) by radiative 
cooling, the temperatures suddenly rise with the shock arrival. Due 
to fast cooling from atomic hydrogen, the post-shock layer cannot 
be resolved until the temperature falls back to ~ 10 4 K, where 
the post-shock pressure balances with the ram pressure. However, 
our insufficient resolution there does not affect the evolution there- 
after because the time-scale for the atomic cooling is very small 
compared with the entire evolutionary time-scale of the shell. For 
element (4) in Fig. [3] the radiative cooling rate and the chemical 
fractions are shown in Figures|4]and[5] respectively. Following the 
shock heating, the gas cools to somewhat below 10 4 K through 
the Ly a emission. By this time, abundant H2 has formed as a re- 
sult of a large amount of available electrons. The H2 abundance 
reaches as high as ~ 10 -3 (Fig. [5j, about an order of magnitude 
higher than that formed in an almost neutral gas. Below ~ 8000 
K, the H2 becomes the most important coolant (Fig. |4j. The H2 
cooling reduces the temperature down to ~ 200 K, giving rise to 
compression of the gas. The metal line cooling has a negligible ef- 
fect during this phase because of the low abundance of metals (Fig. 
|4). Below ~ 200 K, the H2 cooling rate significantly drops and 
the other coolants cannot compensate for the inefficiency. At this 
moment, radiative cooling cannot reduce the shell pressure rapidly, 
while the ram pressure decreases faster than the shell pressure. As 
a result, the shell pressure becomes higher than the ram pressure 
and the shell begins to gradually distend to maintain the balance of 
pressure. This is seen in Fig.[3]as the turn-around of the trajectories 
at ~ 200 K. The adiabatic expansion phase continues for ~ 10 6 - 
10 7 yr, which is much greater than a typical time-scale for the H2 
cooling phase (~ 10 5 - 10 6 yr), until the shell dissolves into the 
ambient medium. 

At low temperatures (< 150 K), HD is known to form abun- 
dantly via the rightward reaction of 

D + + H 2 ^ H+ + HD. (6) 

Although the temperature falls below the threshold for HD forma- 
tion in our case, the density within the shell decreases as well. The 
reaction thus proceeds slower than the evolutionary time-scale of 
the shell, and the HD fraction is frozen to the final abundance of 
~ 10 -6 , which is an order of magnitude smaller than the total deu- 
terium abundance (Fig.[5]l. The gradual density decrease causes the 
HD cooling to be inefficient. For T < Tcmb , the CMB heating via 
HD lines is weak and the trajectories in Fig.|3]exhibit the relation 
for the adiabatic expansion, T oc n 2 ^ 3 , except the highest density 
case of trajectory 1. 

Note that the cooling time in the post-shock gas lengthens for 
smaller shock velocity, and reaches ~ 10 B yr for a gas entering the 
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Figure 1. Evolution of radial distributions of (a) density, (b) temperature, (c) pressure, and (d) velocity for the fiducial run (EgN = 10 erg, no = 1 cm , 
Z = 1O~ 4 Z ). Shown are the distributions at (1) 9.0 X 10 3 , (2) 4.8 X 10 4 , (3) 1.7 X 10 5 , (4) 4.9 X 10 5 , (5) 1.1 X 10 6 , (6) 5.0 X 10 6 , and (7) 1.0 X 10 7 
yr, respectively. The epochs (1) and (2) correspond to the Sedov-Taylor phase, and those from (3) to (7) to the snowplough phase. 



shell during the MCS phase. Thus, the temperature at the outer edge 
of the shell remains high (Fig. [2]b). We use this result to explore 
the conditions for the shell fragmentation in Section [Jl4l 



3.3 Metallicity Effects on the Shell Evolution 

For a high metallicity gas, metal cooling or heating becomes im- 
portant at low temperatures (< 100 K). In this section, we elab- 
orate on this effect on the evolution of the shell. Fig. [6] shows 
the thermal evolution of some fluid elements for the parameters 
Esn = 10 52 ergs, n = 1 cm" 3 , and Z = 10~ 2 Z©, which is 100 
times higher in metallicity than in the fiducial case. The trajecto- 



ries are almost the same as those in the fiducial run until they reach 
~ 200 K where the H2 cooling rate sharply decreases. In this case, 
the Ol and then Cl become the dominant coolants for < 300 K 
and < 100K, respectively (Fig. [7}. Unlike in the fiducial case, the 
temperatures approach the CMB value through radiative cooling 
rather than through the adiabatic expansion. With a rapid decrease 
in the radiative cooling rate near the CMB temperature, the shell 
begins to distend in a manner similar to the fiducial case. Owing 
to the metal cooling, the temperatures at which the shell starts to 
swell are lower than those in the fiducial case. When the tempera- 
ture falls below the CMB value, the radiative heating via Cl lines 
is effective enough to keep the temperature nearly constant at ~ 40 
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Figure 2. Contour maps for densities (a) and temperatures (b) in the shell as a function of time for the fiducial run (Egjj = 10 52 erg, no = 1 cm -3 , 
Z = 1O _4 Z ). The longitudinal axis is the column density measured from the centre. The two thickest lines represent the boundaries of the shell. In 
figure (a), we show the iso-density surfaces of 2000 cm~ 3 (dashed), 1000 cm~ 3 (dotted-and-dashed), 500 cm~ 3 (double-dotted-and-dashed), 100 cm -3 
(dotted), 50 cm -3 (solid), and 10 cm -3 (thick-dashed line). Figure (b) represents the iso-temperature surfaces of 1000 K (solid), 500 K (dashed), 200 K 
(dotted-and-dashed), 100 K (long-dashed), 50 K (dotted), and 25 K (double-dotted-and-dashed line) . 
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Figure 3. The evolution of four fluid elements on the density-temperature 
diagram in the fiducial ran. The initial locations of the fluid elements are 
(1) 68, (2) 77, (3) 88, and (4) 104 pc from the centre. Those elements en- 
ter the shell in the snowplough phase. The arrows indicate the direction of 
evolution. The circles, triangles, squares and black triangles are indicate the 
times 10 5 yr, 5 X 10 5 yr, 10 6 yr and 5 X 106 yr, respectively, after which 
the shock-compressed gas cools to 10 4 K. 



K, although with a gradual decrease in the temperatures. Thus, the 
final temperatures are somewhat higher than in the fiducial case. 

How much metallicity is needed to affect the evolution of the 
shell? To answer this question, we performed other simulations 
with metallicity Z = 10 -3 Zq, and found that the metal cooling 
remains below that of HD molecules. Hence, in this case, the evo- 
lution of the shell is almost the same as in the fiducial case. We con- 
clude that the minimum metallicity necessary to affect the thermal 
evolution of the shell is ~ 1Q~ 2 Zq. We will mention its influence 
on the conditions for the shell fragmentation in Section [3~4l 

3.4 Fragmentation of the Shell 

3.4.1 Overview of Linear Perturbation Analyses 

Here, we discuss whether the supernova-driven shell fragments into 
pieces before it dissolves into the ambient medium. When the shell 
fragments by gravitational instability, the next-generation stars are 
expected to form from the fragments. Conditions for fragmen- 
tation of an expanding and decelerating shell have been studied 
by several authors both analytically (Vishniac 1983; Nishi 1992; 
Elmegreen 1994; Whitworth et al. 1994) and numerically (Yoshida 
& Habe 1992; Mac Low & Norman 1993). The linear analysis 
by Elmegreen (1994) showed that the instantaneous growth rate 
of density perturbations for transverse motion at wavenumber k is 
given by 



V V 2 
«(M) = -3 5 + VlP 



fc 2 c 2 + 27rGE fc, 



(7) 



where R, V, En, and c s are the radius, velocity, unperturbed column 
density, and sound speed of the shell, respectively. Here, the per- 
turbations are assumed to grow in a way similar to exp(u)t). The 
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Figure 4. Radiative cooling and heating rate of a typical fluid element (4 in 
Figure [5) against the temperature. The vertical dotted line shows the CMB 
temperature, T^mb = 44 K, below which the processes work as heating. 
The lines represent the radiative rates of H2 (long-dashed), HD (dashed), 
Cl (dotted) and Oi (dotted-and-dashed). The solid line shows the total rate, 
including the H atomic cooling. 
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Figure 5. Evolution of the H2 (solid), HD (dotted-and-dashed) and electron 
(dotted) fractions at the same fluid element as Figure [4] against the temper- 
ature. The horizontal dashed-line represents the total deuterium abundance 
of 2.5 X 10~ 5 . The effect of the shock is observed as the kinks in the frac- 
tion curves. Thereafter the temperature continues to decrease, while the H2, 
and HD fractions gradually increase. 
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Figure 6. Same as Fig. [5] except for the metallicity Z = 10 2 Zq. The 
horizontal dashed-line represents the CMB temperature at z = 15. 
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Figure 7. Same as Fig. [4] but for the fluid element (4) in Figure [6] with 
metallicity 1O~ 2 Z . 



seeds of density perturbations are generated, for instance, by the 
decelerating shock instability (Mac Low & Norman 1993), or by 
the shock propagation into the inhomogeneous medium. The per- 
turbations grow by collecting the ambient medium, while the new 
mass accretion onto the shell increases total transverse momentum. 
In addition, the spherical expansion stretches and attenuates pertur- 
bations. These processes are taken into account via the V/R terms 
in equation ifTJ and the net effect acts to hinder growth of the per- 
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turbations. The maximum growth rate, 



+ 



ttGEq 



is attained at wavenumber 



ttGEq 



The shell becomes gravitationally unstable if oj max > 0, i.e.. 



te 



R ^ 2 c s 
V > ttGEq 



+ 2 

i a 



tc 



(8) 



(9) 



(10) 



where t cxp , tff, and t C ro SS are the shell expansion time, free-fall 
time in the shell, and the sound crossing time for the shell width, 
respectively. 

Even when the shell becomes gravitationally unstable, it takes 
some time for the perturbations to grow to induce the fragmenta- 
tion. The time-scale for fragmentation can be roughly estimated 
by l/w m ax. It should be noted, however, that the growth rate for a 
fixed wavenumber changes with time because V/R and Eo in equa- 
tion are modified by the mass accretion. To evaluate the typical 
fragmentation scale, we keep track of all the modes that attain the 
maximum growth rate during the shell history. For those modes, we 
define the fragmentation time tf as: 



oj(k,t)dt = 1, 



(11) 



where ti is the time when the mode k first becomes unstable. 
Namely, tf is the time for the perturbation to grow e times the 
initial value. When condition dl It is met, the shell is assumed to 
fragment into pieces with scales corresponding to fc max (Ehlerova 
et al. 1997; Elmegreen, Palous & Ehlerova 2002). We evaluate the 
fragment mass by 

M frag (A max ,t / )=7rf^ N ) Eo(t/), (12) 



where A max = 27r/fc max . 

The shell loses its identity when it slows down to the sound 
speed of the ambient medium. For triggering subsequent star for- 
mation, fragmentation must occur within the lifetime of the shell. 

3.4.2 Application of the Linear Analysis 

In this section, we use the result of the linear analysis to find the 
conditions for fragmentation of a supernova-driven shell. The dis- 
persion relation l[8} is derived with one-zone approximation to the 
shell, whereas the shell has a stratified structure in our simulations. 
The application of the linear analysis to our result is, therefore, not 
straightforward. If the stratified structure is maintained until the 
shell disappears, the gravitational instability seems to occur first 
in a cold layer of the shell due to its low temperature. However, 
this is not the case because the mass of the cooled gas is not suf- 
ficient to trigger the gravitational instability. On the contrary, the 
coldest layer tends to gradually expand as discussed in Section [X2l 
indicating that the self-gravity is not important there. In reality, the 
gravitational instability most easily occurs at the scale comparable 
to the width of the shell because the column density of the layer 
must be high enough for fragmentation. We then simply estimate a 
representative temperature by taking mass-weighted average over 
the whole shell. 

In Fig. [8] we show the parameter range of (-Esn, no) where 
the shell fragments within its lifetime for runs with Z = 10 -4 Zq 



and 10 -2 ' Zq. Metallicity affects the result only for (Esn, no) = 
(10 52 , 10), where the 1Q~ 4 Zq run satisfies condition j 10b but not 
dl lb before the dissolution of the shell. For high ambient density, 
the shell accumulates a large amount of gas within the lifetime and 
fragments even in the case of W~ 4 Zq. Thus, the metallicity de- 
pendence disappears for such high ambient densities. In the cases 
of Es~n = 10 52 erg, the shell fragmentation occurs for both metal- 
licities at the ambient density ^ 10 2 cm" 3 . For lower explosion en- 
ergy Es~n = 10 51 erg, this threshold density increases to 10 3 cm -3 . 
In Fig. [8] we also plot the conditions for shell fragmentation found 
by Salvaterra et al. (2004) and Machida et al. (2005) for compar- 
ison. The discrepancy among their conditions stems from the dif- 
ference in their adopted shell temperatures and criteria for frag- 
mentation. Our result almost coincides with that by Salvaterra et 
al. (2004) because their assumed temperature is almost the same as 
ours. As long as the mean temperature falls to 200 - 300 K through 
H2 cooling, the total amount of gas swept up into the shell does 
matter for fragmentation. Thus, the conditions for fragmentation 
rely heavily on the ambient density as well as explosion energy, but 
little on the metallicity. 

In Table Q] we show the properties of the shell at fragmenta- 
tion time tf. There is a range of the fragment mass because some 
unstable modes fulfill the condition dl lb soon after t / . This means 
that the fragmentation scale depends largely on the shape of initial 
perturbation, or that the shell fragments with a wide range of mass 
scales. The fragment mass for models with no 10 2 cm -3 is in the 
range 10 2 - 10 3 Mq, while that for a model with no = 10 cm -3 is 
an order of magnitude lower. Because the wavelength A max length- 
ens for the smaller column density (see equation^, the fragment 
mass becomes higher for lower ambient density. 



4 DISCUSSIONS 

4.1 Significance of Star Formation Triggered by Supernova 
Explosions 

The supernova explosion is an important feedback process influ- 
encing cosmic star formation history as well as the nearby star for- 
mation activity. We have shown that supernova explosions induce 
subsequent star formation as a result of shell fragmentation in a 
high-density environment (see Fig. [8}. Here, we estimate the criti- 
cal mass of haloes above which second-generation stars form in the 
induced manner, using our results. We assume that the explosion 
occurs at the centre of the halo. We have derived the minimum am- 
bient density p m i n for a given explosion energy, at which the shell 
fragments at a radius n. Thus, the condition for the shell fragmen- 
tation is that the column density exceeds p m i n rf /3 at a radius < rt , 
or in other words, the mean density pi within the radius n is lower 
than p m i B . We evaluate the mean density pi assuming the density 
profile in the halo p(r) cx r~ w , with the index w ~ 2 as suggested 
by numerical calculations of the first star formation. 

Using the mean density p v i I (z c ) and virial radius r v ir of a halo 
viliarizing at z c , we obtain the mean baryonic density inside rf as 



Pc 



O-Pvir(Zc) ( — 
a &771 \ / vir 



3.0 x 10" 25 



gem 



l + z c 
16 



Tvit 



(13) 



(14) 



Here, we apply the spherical collapse model to the mean density in 
the halo, p v \ r — 18TY 2 (Qt/Q,rn)pb(z c ), where pt, is the background 
density (e.g., Padmanabhan 1993). By eliminating r v i r in equation 
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Table 1. Properties of the shell when the shell fragmentation first occurs. The left column is the model parameters. The notation T SI1 represents the mass- 
weighted mean temperature in the shell, So the column density, AR the width of the shell, p sh the mean density in the shell, t max the time when the shell 
becomes gravitationally unstable first, 1 1 the fragmentation time, A max the wavelength of perturbation which becomes unstable first, i? s i n and M sn the radius 
and mass of the shell at fragmentation, and Mf rag the mass of the fragments, M the minimum mass of the fragments after tf , the maximum mass 

of the fragments after tf . The last column is the critical halo mass below which star formation cannot occur in the induced manner. 
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Figure 8. Dependence of the fragmentation result of the shell on the am- 
bient density and explosion energy for metallicities Z = 10 — 2 Zq and 
Z = 10~ 4 Z Q . In the cases of E SN = 10 51 and 10 52 erg, the circles, 
triangle, and crosses mean that the fragmentation occurs for both metallic- 
ities, only for Z = 10~ 2 Zq, and for neither, respectively. In the cases of 
^SN = 10 53 erg, the circles and crosses represent the successes and fails 
of the fragmentaion respectively, and the triangle satisfies the condition )10t 
but does not equation {TT}. Also plotted are the results by Salvaterra et al. 
(2004; solid) and Machida et al. (2005; dotted), above which the fragmen- 
tation occurs. 



d!3t with A/haio = (47r/3)pvi r rvi r , the halo mass Mhaio can be 
expressed as 



■ * 47T _-172, s 3 I £lb 



-3/2 



-3/2 
Pi ■ 



(15) 



1.5xlO B M 







l + z c 
16 



n 



\10 2 cm- 3 / 



Setting pi = p m in, we obtain the critical halo mass for triggered 



41.4pcy Vl0 2 cm- 3 / 
(16) 

where n m m = pmin/ ftrriu- Note that the dependence on ambient 
density n m i n is also included through the fragmentation radius rf. 
The critical halo masses range from 1.8 x 10 7 M to 9.2 x 10 s M Q . 
These values are listed in TableQ] These mass-scales correspond to 
ss 3cr fluctuations virializing at z = 15 and the abundance of such 
objects is about a few per cubic Mpc. The critical mass for trig- 
gered star formation M{^i exceeds the halo mass of 4.7 x 10 7 
M [(1 + z c )/16] _3//2 corresponding to the virial temperature of 
10 4 K(Barkana & Loeb 2001), and consequently, star formation 
triggered by a supernova explosion is unlikely to occur in low- 
mass haloes (T v i r < 10 4 K) where primordial stars form solely via 
H2 cooling. Instead, in such low-mass haloes, the expanding SNR 
evacuates all the gas and quenches the subsequent star formation 
activity. The triggering process becomes increasingly important at 
lower redshift as more massive haloes of > 10 8_9 M Q begin to 
form. 

So far, we have assumed the density distribution p ot r~ 2 , 
which is appropriate at the stellar birth. If an expanding Hll region 
created during the stellar lifetime disturbs and even blows the sur- 
rounding material away from the host halo before the supernova 
explosion, the column density of the supernova-driven shell is sig- 
nificantly reduced and the subsequent triggering process does not 
work. This effect is important in low-mass (< 1O 6 M ) haloes (Ki- 
tayama & Yoshida 2005; Whalen et al. 2008). However, in massive 
haloes (> 10 s M©) where the triggered process is possible, the Hll 
region is trapped deep inside the virial radius and barely removes 
any gas from the host haloes. Thus, the expansion of the Hll region 
does not change our criteria for the triggered star formation. Here, 
it should be noted that the densities in the centre of the trapped 
Hll regions are 10 4 - 10 7 cm -3 and much higher than we have 
examined, according to some numerical simulations (Kitayama & 
Yoshida 2005; Whalen et al. 2008). The reason why we have inves- 
tigated the density range of 10 2 - 10 3 cm 3 is only for comparison 
with the previous works of Salvaterra et al. (2004) and Machida et 
al. (2005). 

Whalen et al. (2008) studied supernova explosions in neutral 
haloes undisturbed by the expansion of Hll regions. They showed 
that the SNR is trapped by the deep gravitational potential well 
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of the dark matter. It initially expands out to ~ 50 pc, and then 
turns to fall back. Metals in the ejecta will mix with gas in the 
halo before the SNR stagnates. The polluted gas may fragment 
into clumps through metal fine-structure line cooling, which leads 
to subsequent star formation. This triggering mode is beyond the 
scope of this paper, and to be examined in future studies. In addi- 
tion, we have ignored another process of possibly triggering star 
formation described below because of a limitation of our simula- 
tions. The Hll region prior to the supernova sweeps up materials 
into a shell (Kitayama et al. 2004; Whalen et al. 2004). The SNR 
catches up with the pre-existing shell. Whalen et al. (2008) showed 
that the collision is quite violent and would mix heavy elements 
from the ejecta with materials in the shell. Fracturing the ejecta 
into the dense shell and the resultant enhanced cooling could lead 
to star formation. The shell/ejecta interaction occurs even in low- 
mass halos of < 1O 8 M0. Thus, simulations with more realistic 
density distribution are needed for a complete survey of the condi- 
tions required for triggered star formation. 

We finally note that hydrodynamical mixing of metals in the 
SNR is generally important for realistic evolution. Even if the am- 
bient medium is constituted by pristine H and He gas, the polluted 
gas will cool via metal fine-structure lines. Although our calcu- 
lations with pre-enriched ambient medium predict that the criti- 
cal halo mass needed for triggering star formation depends hardly 
on metallicity, it is still uncertain how the hydrodynamical mixing 
spreads metals into the SNR and affects the triggering processes. 

4.2 Subsequent Evolution of the Fragments and Possibility of 
Low-mass Star Formation 

After the fragmentation, each fragment continues to contract 
through its own gravity. Although the mass of fragments is ~ 10 3 
- 10 4 Mq (TablefTJ, they are expected to fragment once again into 
smaller (< 10 Mq) pieces. The shell initially fragments into disk- 
like clouds, and then further fragmentation leads to smaller fila- 
mentary clouds (Miyama, Narita & Hayashi 1987). Nakamura & 
Umemura (2001, 2002) examined fragmentation of primordial fil- 
amentary clouds and concluded that they fragment into spherical 
cores of 1 - 10 2 Mq, the exact value of which depends on the ini- 
tial cloud density n c and gravity-to-pressure ratio /. The fragmen- 
tation leads to low-mass (~ IMq) cores only if the conditions that 
n c > 10 5 cm" 3 and / > 3 are met. Salvaterra et al. (2004) pointed 
out that these conditions are easily satisfied and thus low-mass star 
formation occurs in the supernova-driven shell in primordial envi- 
ronments. However, our more detailed calculations demonstrated 
that the number density in the shell never exceeds 10 5 cm" 3 for 
the ambient density no ^ 10 3 cm" 3 , and the low-mass stars are 
unlikely to form without any metal enrichment even in the triggered 
fashion. 

Provided that the metallicity is higher than Z CI i t = 10" J - 
10~ 6 Zq and a sizeable portion ~ 0.5 — 1 of those metals condenses 
into dust grains, fragmentation is induced by dust cooling at high 
densities > 10 10 cm" 3 (Schneider et al. 2002, 2006; Omukai et 
al. 2005; Tsuribe & Omukai 2006, 2008; Clark, Glover & Klessen 
2008). In the early universe, dust grains originate from supernovae 
(Todini & Ferrara 2002; Nozawa et al. 2003; Schneider, Ferrara & 
Salvaterra 2004) rather than the AGB stars. Nozawa et al. (2007) 
studied the formation of dust grains in population III supernovae 
and their destruction by the reverse shock. They concluded that the 
mass of the surviving dust grains in the shell is at most 0.1 Mq for 
the ambient density no = 10 cm -3 at t ~ 0.7 Myr. To reach the 
critical metallicity Z CT , the total mass in the shell should be less 



than 5 x 10 6 (Z cr /10~ 6 Zq)~ x M before fragmentation, under 
the assumption of homogeneous grain distribution within the shell. 
Because the shell mass at fragmentation is ~ 10 6 Mq for no = 10 
cm" 3 (Sec[3}, low-mass stars can be triggered by the first super- 
nova. 

4.3 Effects of External Photo-dissociating Radiation 

We assume that H2 molecules in the ambient medium are com- 
pletely photo-dissociated by the central star before the supernova, 
but they are replenished quickly due to a high fraction of the cat- 
alyst electrons in our calculations. As a result, the ambient tem- 
perature decreases to several hundreds of Kelvins owing to the H2 
cooling (Fig. []}. However, we have not considered the effects of 
any external radiation. Once stars appear in the universe, they build 
up a background of H2 dissociating radiation, which delays subse- 
quent star formation (e.g., Haiman, Rees & Loeb 1997; Haiman, 
Abel & Rees 2000). If an external radiation flux dissociates the 
H2 molecules, the gas temperature remains higher because of inef- 
ficient cooling. The low temperature in the ambient medium pro- 
longs the lifetime of the shell. The fragmentation time is slightly 
shorter than the shell lifetime even in our calculation, and thus the 
fragmentation may not occur with external heating. 

According to a cosmological simulation by Johnson, Greif & 
Bromm (2008), the photo-dissociating radiation flux is J21 < 0.04 
at 12.87 eV in units of 10" 21 erg s" 1 cm" 2 Hz" 1 sr" 1 . For 
this background, the H2 dissociation timescale is tdiss ~ 6 x 
10 5 ( J21/O.O4)- 1 yr (Abel et al. 1997), whereas the typical H 2 for- 
mation time behind the shock in our calculation is t{ olm ~ 10 4 
yr at the beginning of the PDS phase. The formation time is thus 
an order of magnitude shorter than the dissociation time. Hence, 
the external photo-dissociating flux does not prevent the H2 for- 
mation at this phase. Soon after that, the H2 column density of 
the shell rapidly increases to shield the photo-dissociation (i.e., 
> 10 14 cm -2 ; e.g., Draine & Beltoldi 1996) as seen below. Before 
the photo-dissociation, the H2 column density of the shell reaches 
~ y(H2)nHf s h£diss, where ?/(H2) and v s h are the fractional abun- 
dance of H2 molecules and shock velocity. The shock velocity is 
on the order of 10 2 km s" 1 at the incipient PDS phase, while hh 
is enhanced to ~ M 2 no by an isothermal shock with the Mach 
number A4. The H2 fract ion necessary for self-shielding is then 
only 

in - 8 ( _^h_v 1 ( tdiss y 1 ( n ° y 

vh 2 ~iu y W 2 km/s J \ W Gyr) VI cm" 3 / ' 

(17) 

where we assume that M ~ 10. Because the H2 abundance far ex- 
ceeds this value within the dissociation time (Fig. [5}, we conclude 
that the external photo-dissociating radiation does not significantly 
prevent H2 formation in the shell. 

4.4 Effects of Non-gravitational Instabilities 

In our analysis of gravitational instability, we have implicitly as- 
sumed that the coherence of the layer is not destroyed by some dy- 
namical instabilities emerging before the gravitational one. For in- 
stance, the decelerating shock instability is expected to grow from 
the linear analysis (Vishniac 1983). Mac Low & Norman (1993) 
numerically demonstrated, however, that this instability saturates 
at a rather low amplitude in the non-linear stage, and does not sig- 
nificantly disturb the shell structure. The thermal instability is also 
expected to occur in the rapid cooling phase in the shocked layer. 
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Koyama & Inutsuka (2002) studied propagation of a shock into 
the contemporary ISM and fragmentation of a shock-compressed 
layer, taking account of cooling, heating, and thermal conduction. 
They showed that the thermal instability causes the cooled layer to 
fragment into small (~ 10 - 100 AU) cloudlets with some trans- 
lational motions. Bromm, Yoshida & Hernquist (2003) simulated 
the first supernova explosions, and showed that the swept-up shell 
consists of small cloudlets, which are presumably formed by ther- 
mal instability. The motions of the cloudlets effectively enhance 
the sound speed within the shell, leading to delay in growth of the 
gravitational instability. After the shell fragments through the grav- 
itational instability, the cloudlets will aggregate into a larger cloud 
with the contraction of the fragments. This may ultimately lead to 
cluster star formation (Elmegreen 1998). Thus, thermal instability 
influences not only the shell fragmentation but also the subsequent 
star formation in the fragments. Although our one-dimensional cal- 
culations cannot clarify the detailed fragmentation processes in 
the shell and subsequent evolution of the fragments, future multi- 
dimensional simulations are expected to reveal them. 



5 SUMMARY 

To explore the possibility of star formation triggered by super- 
novae in the early universe, we have studied the evolution of shells 
formed around the supernova bubble during the snowplough phase 
in low-metallicity environments. We investigated detailed structure 
of the shell, using spherical symmetric hydrodynamics with non- 
equilibrium chemistry for different sets of ambient densities, ex- 
plosion energies, and metallicities. Our results are summarized as 
follows: 

(i) We have demarcated the ranges of metallicity in terms of ther- 
mal evolution of the shell: 

- For metallicities with < W~ 2 Zq, the post-shock gas first 
cools promptly to ~ 200 K through the H2 cooling, and then 
progressively to a lower value through adiabatic expansion. Al- 
though 10 % of the deuterium transforms into HD in the shell, 
the HD cooling and heating have never become important be- 
cause of decreasing density. 

- For higher metallicities (> 10~ 2 Zq), the fine-structure line 
transitions of O and C reduce the post-shock gas below 100 K. 
After the radiative cooling becomes inefficient, the gas then adi- 
abatically expands as in the lower-metallicity cases. The CMB 
heating via the Cl fine-structure transitions hinders the gas from 
cooling below the CMB temperature. 

(ii) We have examined whether the swept-up shell becomes grav- 
itationally unstable before it mixes with the ISM, using the linear 
perturbation analysis. Fragmentation of the shell easily occurs in 
denser ISM and for higher explosion energies because these con- 
ditions are favourable to attaining the column density of the shell 
necessary for fragmentation. The criteria for the fragmentation de- 
pend little on metallicity as thermal evolution of the shell is not so 
sensitive to it. 

(iii) We have evaluated the fragment mass from the wavelength of 
the maximum growth modes. Because the wavelength is inversely 
proportional to the column density, the fragment mass becomes 
lower for both higher ambient densities and higher explosion en- 
ergies. The resulting fragment mass is 10 2 - 10 3 M Q . 

(iv) We have derived the critical halo mass M"jj ~ 10 s M 
above which star formation can be triggered by a supernova explo- 



sion. In minihaloes (< 10 Mq) where the first stars are more likely 
to form, this triggering mechanism does not work. 
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APPENDIX A: CHEMICAL REACTIONS 

We list the chemical reactions we include in Table I All 

The definitive version is available 
' www.blackwellsynergy.com' . 



Star Formation Triggered by Supernova Explosions in Young Galaxies 



Number Reactions References 
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Table Al. References.— 1. Abel et al. (1997); 2. Galli & Palla (1998); 3. 
Shapiro & Rang (1987); 4. Stancil, Lepp, & Dalgarno (1998); 5. Galli & 
Palla (2002); 6. Palla, Salpeter & Stahler (1983); 7. Stancil et al. (1999); 8. 
Nahar & Pradhan (1997). 



